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Abstract 

By introducing the divergence of a vector potential into the Lagrangian, a Lagrangian framework is 
developed to incorporate the surface energy into elasticity. Besides the Euler-Lagrange equation and 
natural boundary condition, a new boundary constitutive equation is derived from the variation of the 
Lagrangian and configuration on which the Lagrangian is defined. On the boundary surface, explicit 
expression of the vector potential with respect to field variable and surface curvature is determined. 
Based on this framework, an elastic model with relevance to the surface energy is established. The 
Young-Laplace's formula is generalized into elastic solid in a new form. Making use of this model, 
we investigate the surface energy effect in the radial vibration of spherical nanoparticle. Numerical 
calculation shows that natural frequencies of nanoparticle will shift down due to the surface energy. 
This shift is especially apparent in the vibration of soft matter nanoparticles. 

Key words: Euler-Lagrange equation, surface energy, generalize Young-Laplace' formula, natural 
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1 Introduction 

From the view of atomic level, the free surface of a bulk material is different from its interior. This 
difference gives rise to a excess energy to appear on surface. It is called the surface free energy (or 
shortly called the surface energy). In classical continuum mechanics, the influence of surface energy on 
the deformation of solid is usually ignored, because atoms on surface is far smaller in magnitude than 
that in the interior of bulk. However, with entering from microscale into nanoscale, the ratio of atoms on 
surface to total atoms becomes gradually significant, therefore, the surface energy effect has to be taken 
into account in the deformation of solid. 

On the surface energy effect in elastic deformation, early works can be traced back to Shuttleworth 
[1], Herring [2], Nicolson [3] and Vermaak et al [4]. They investigated the lattice contraction caused 
by the surface tension in small crystal particles. A general theoretical framework incorporating surface 
energy into continuum mechanics was proposed in the 1970s by Gurtin and Murdoch [5, 6]. Within 
this framework, Steigmann and Ogden [7] generalized the conception of surface energy to the curvature 
dependence case. Huang and Wang et al [8, 9] clarified the necessity to introduce three configurations in 
an accurate and complete description for the deformation of continuum with surface energy. Altenbach 
et al [10] gave a proof on existence and uniqueness of the solutions of initial and boundary value problems 
of linear elasticity with surface energy. Schiavone and Ru [11] discussed the solvability of boundary value 
problems in a theory of plane-strain elasticity concerned with surface energy. 

Various models based on the Gurtin-Murdoch theory have been presented in several contexts, e.g. 
nanoparticles, wires and films (Huang and Thomson et al [12]; Cammarata [12, 13]; Dingreville and 
Qu et al [15]; Park [16]), nanoscale rods, beams, plates and shells (Miller and Shenoy [17]; Huang [18]; 
Bar and Altus et al [19]; Liu and Rajapakse [20]; Wang and Feng [21]; Wang and Zhao et al [22]), 
nanoinclusions (Duan and Wang et al [23]; He and Li [24]; Sharma and Ganti [25]). So far, a considerable 
number of literatures have been devoted to the continuum models considering surface energy. It is not 
the purpose of this paper to list and review these abundant literatures. The reader can be referred to the 
recent reviews by Wang and Huang et al [26] and Duan, Wang and Karihaloo [27] on the literatures. 



The subject of the paper will focus on developing a Lagrangian formulation corresponding to the Gurtin- 
Murdoch theory. The novel feature of this formulation consists in that a vector potential is introduced 
to characterize the surface energy of continuum. Consequently, the Eular-Lagrange equation, natural 
boundary condition and a new boundary constitutive equation can be consistently determined by simul- 
taneously taking the variation of the Lagrangian and configuration of continuum on which the Lagrangian 
is defined. 

The paper is outlined as follows. In section 2, based on a vector potential, a Lagrangian framework is 
developed to incorporate surface energy into continuum mechanics. The boundary constitutive equation 
and the vector potential are determined on the boundary surface. By linearization, in section 3 we put 
forward to an clastic model with relevance to surface energy in which a generalized Young-Laplace's 
formula is derived consistently. Making use of this elastic model, we investigate the surface energy effect 
in the radial vibration of spherical nanoparticle. some interesting results are predicted. Finally, we close 
this paper with a summary. 

Notation: A compact notation is used, with boldface letters being vectors or tensors. The index rules 
and summation convention are adopted. Latin indices have the range 1, 2, 3. Partial derivatives with 
respect to coordinates are represented as dk = djdx^ or (•)^fc = d{-)/dx''. Partial derivative with respect 
to time is denoted by an upper dot, e.g., a = da/dt. Other symbols will be introduced in the text where 
they appear for the first time. 

2 Theoretical framework 

Let X = {x^} (j = Ij 2, 3) be a 3-dimensional position vector on C i?^. The Greece letter J7 stands for 
current configuration of a continuum. We use (f = (fi{t,x.) to characterize the field defined on [to,^i] U fl, 
where t G [to,ti] denotes time. Depending on circumstances, <^ is a scalar, vector or tensor. 
As is well known, if two Lagrangians differ by the divergence of a vector function only of the field variables 
and coordinates, their Euler-Lagrange equations are the same in form. Thus, if such a divergence term 
is added to a Lagrangian, dynamic behaviors of field are unchanged under the fixed boundary condition. 
However, if the boundary is free, the additive divergence will contribute to the field by entering the 
natural boundary condition. Therefore, a problem worthy of being asked is whether a divergence term 
added to the Lagrangian can be used to characterize the surface effects of a free continuum or not. 
Motivated by this idea, we assume a Lagrangian can be written as L = L{(p, (p,dj(p) — dkY''{x, (p), where 
yk — yfc(^x, if) denotes a vector potential. As thus, the action functional of tp can be represented as 

A[<p] = r f [L{<p, cp, djip) - dkY'^i^, ip)]dv{^)dt. (1) 

Jto Jn 

In terms of the divergence theorem, A[<p\ can be also written as 

A[ip] = r j L{^,cp,dj<p)dv{^)dt- r [ F'=(x,<pKda(x)di, (2) 

where nk denotes the unit normal vector on the boundary dil of f7. It is necessary to emphasize the dis- 
tinction between the vector potential and surface energy. The former is defined on continuum, including 
its interior and surface; but the latter is only on the surface. In the following, we will demonstrate that 
they are equivalent on surface. 

If surface energy is concerned, the changes of the configuration of continuum will become a factor having 
to be considered. Consequently, the calculation for the variation of A[lp\ shoiild not only include the 
variation of the Lagrangian but also the variation of the surface area and volume of continuum. Under 



this consideration, the variation of A[ip] can be calculated as follows: 
SA[ip] = [ f (5Ldw(x)dt+ / / Ldv{x.)dt 

Jto Jn Jto J so, 

- r I 5Y^nkda{^)dt ~ T I Y^nkda{^)dt 

Jta Jdn Jto Jsian) 

- / -jr—nkSipda{x)dt + / / Ldw(x)di - / / Y''nkda{x)dt 
Jto Jan (^V Jto Jsn Jto Js{an) 

f*' f ,dL d,dL^ 9L , , f ,dL 



n f ,dL d.dL. ^ , dL ^ , /"r^^c It, J . X 

f . dL dY\ , J , 

^Jto iaat^ - 



/to 

+ [ ' [ Ldv{yi)dt - [ ' f y'=nfcda(x)dt, (3) 
^to ^50 .^to Js{an) 



in which we have supposed that <p{x), L{(p, (p, dj(p) and y (x, are suitably smooth functions. Ou-Yang 
and Helfrich [28] have proved that 

/ Ldv{x) = I LnjSu^da{x). (4) 
Jsn Jan 

[ y'=nfeda(x)= [ 2HgkjY^5u^da{x), (5) 
J6{dn) Jan 

where 5u^ is the variation of displacement of the area element da(x) on the boundary surface dCl, H 
the mean curvature at x on dVl and gkj the materic tensor. Sip = at the two end points of time, so 
substituting Eq. (4) and (5) into (3) yields 

SA[ip] = / / - ^(^) - d,ij^^)]5ipdv{x)dt+ / / [j^^ - ^]nkS^da{x)dt 

+ 11 {Luj - 2HgkjY^)5u^da{x)dt. (6) 
Jta Jan 



According to the Hamilton's principle, 5A[ip] = 0. Therefore, the fundamental lemma of variation leads 

to the below results: 
Euler-Lagrange equation: 

d ,dL . ^ , dL ^ dL ^ 
dM^^'^^didwi'-d^^'- 



Natural boundary condition: 

, dL dY'', 

\nk 



= 0. (8) 

90 



d{dk(p) d(p 
Boundary constitutive equation : 

{Ln,-2HgkjY'')\g^^ = 0. (9) 

Eq. (7) and (8) show that a divergence term, only relevant to the field variable, added to the Lagrangian 
has no effects on the Euler-Lagrangc; equation, but it contributes to the natural boundary condition of 
field. Since Y'^ depends only on p, H in Eq. (9) is necessarily relevant to d,y(p. Or else the derivative of 
L with respect to dk^ will be zero, in terms of Eq. (9). That contradicts with Eq. (8). Making use of 
Eq. (9), we have 

, dL „ dH 



-Uj - 2gkjY'' 



= 0. (10) 

an 



'd{d,cp) ' d{d,v) 

Comparing Eq. (10) with (8), one can see that H is linearly dependent of d„ip. Subtracting Eq. (9) from 
(10) gives 

M'' . ..I. dH 



= 0. (11) 
an 



Let denote the contravariant component of the unit normal vector on dCl. So n^nk = S^.. By this 
identity, Eq. (11) can be rewritten as 



= 0, (12) 

on 



where 

BIT 

On dQ, the solution of Eq. (12) can be represented as 

r = 7exp(S<^), (14) 
where 7 is an integral constant. Inserting Eq. (13) in (14) leads to 

y^n.=7exp[2<^n.^^], (15) 

which characterizes the mathematical form of Y'^rik on the free surface of fl. Eq. (15) shows that, on 
surface, Y'^n^ (i-c, F) is a surface curvature-dependent energy. So it is just the surface energy. It should 
be noticed that, due to concerning the curvature of surface, Eq. (15) holds only on the surface 90, but 
not valid in the interior of O. 



3 Linear elastic model with relevance to surface energy 

Consider a free elastic body concerned with surface energy. The Lagrangian of it is assumed to take 
the following form (In the following, all indices are written as subscripts without distinguishing the 
superscripts and subscripts.): 

L = L - Yk^k = -^pukiik - -^CijkiUijUk,i - lfc,fc. (16) 

where Uk denotes the elastic displacement field, p and Cijki are the mass density and elastic tensor. 
Based on the viewpoint of Huang and Wang [8], Yk^k can be regarded as elastic energy transformed from 
surface energy relative to the so-called fictitious stress-free configuration. 
In terms of Eq. (14), YkUk is given by 

YkUk = r = 7exp(S'fcMfc). (17) 
Under the small displacement condition, Eq. (17) is expanded to the second order term of Uk, i.e., 

YkUk = 'y{l + SkUk + ^SkSjUkUj), (18) 

where Sk is determined by Eq. (13)2- Just mentioned above, the mean curvature H in Eq. (13)2 is 
linearly dependent of Uk,i- As thus, H is only represented as H = Ho{l + Uk,k) in terms of the invariance 
of H [29]. Here, Hq is the initial mean curvature. Making use of Eq. (13)2, we have Sk = 2nkHo. With 
help of this result, Eq. (18) becomes 

Yknk = 7(1 + ^HoUkUk + 2HluknkUjnj). (19) 

By means of Eq. (16) and (19), Eq. (7) and (8) can be written as 

CijkiUkji = pui, in f2. (20) 

CijkiUk,inj - 2"fHoni{l + 2HoUknk) = 0, on dfl. (21) 

Eq. (20) and (21) are the equation of motion and boundary condition to characterize the free elastic 
body concerned with surface tension, respectively. By means of the Hooke's law aij = CijkiUk,u Eq. (21) 
is rewritten as 

UijUj = 2'yHoni{l + 2i^o'Wfe^^fe)• (22) 

If the term 2HoUknk is dropped, Eq. (22) will reduce to the Young-Laplace's formula that describes the 
surface pressure of fluid. Therefore, Eq. (22) can be regarded as a generalization of the Young-Laplace's 
formula in solid, whereas 7 represents the specific surface energy factor. 



4 Application to the radial vibration of nanoparticle 



4.1 Chciracteristic frequency equation 

As an example, we consider the radial free vibration of a spherical nanoparticle with the radius of tq. Let 

material be isotropic. The equations governing this vibration can be determined by Eq. (20) and (21) 
together with Hq = — 1/ro and Cijki = \Sij5ki + ^i{SikSji + SuSjk)- In a spherical coordinate system, they 
are written as 

d'^Ur 2 dUr 2 _ 1 d^Ur 

Qj.2 r dr ^ ^2 g^2 

[(A + 2m)^ + 2(A - ^)^Uro = --■ (24) 
or r r 



where A is the Lame constant and fi the shear moduhis. c = ■^/(A + 2/i) / p, being the wave velocity. 
Assume that the solution of Eq. (23) subjected to (24) to have the form below: 

Ur{r,t) = Ur+gir)e"^\ (25) 

where Ur reads 

^'- = -(3A + 2M)ro-V' ^''^ 

which is the radial contraction displacement caused by the surface enrgy [12]. It should be noted that 
Mr > or — >■ 00 whcu rp < Tc = 47/ (3A + 2/i). That is impossible in physics because it indicates that 
the surface enrgy will cause the radial extension of a spherical particle. Therefore, the elastic model is 
no longer valid for the case that ro < Tc- Inserting Eq. (25) in (23) and (24) leads to 

d'^g 2dg 2^ 

[(A + 2Ai)|^ + 2(A-^)^],=.„=0. (28) 
or r r 

Let k = w/c and g{r) = f{r)/y/r. By means of these substitution, Eq. (27) can be transformed into a 
standard Bessel equation. Therefore, the general solution to Eq. (28) can be given as follows: 

g{r) = Ar-'/^J^^2{kr) + Br-'^^Yy^ikr), (29) 

where J3/2{kr) and i3/2(A:r) are first kind and second kind Bessel function of order 3/2, respectively. A 
and B are two undetermined constants. Noticing g{0) = 0, we have B = 0. Thus, Eq. (29) reduces to 

5(r) =^r-i/2j3/2(fcr)- (30) 

Substituting Eq. (30) into (28) yields 

kJ'sMkro) + i^^^^^l^^-^ J3/2(fcro) = 0, (31) 

where J 3/2 denotes the derivative of J3/2 with respect to its argument. Eq. (31) represents the char- 
acteristic frequency equation of the radial vibration of the spherical nanoparticle involving the surface 
energy effect. 



4.2 Results and discussion 

Using Eq. (31), we calculate the natural frequencies of two different nanoparticles. The first is the 
silica gel particle, and the second is the palladium (Pd) particle. The silica gel is a soft matter. We 
take its E = 2.14MPa, 0.48, p 1200kg/m^ and 7 = 22.0mN/m. In terms of the transformation 
A = Eu/[{1 + i/)(l - 2iy)] and = E/[2{1 + 2i/)], we have A = 17.35MPa and /x = 0.72MPa. As thus, 
Tc = 47/(3A + 2/x) = 1.65nm. When calculating, the range of ro should be greater than 1.65nm. 
The changes of first three natural frequencies with the radius ro arc shown in Figure 1. It can be seen 
that the three natural frequencies increase with a decrease in the radius. This is a well known size effect. 
More significantly. Figure 1 illustrates that, with entering the nano-level, the influence of surface energy 
on the natural frequency becomes apparent gradually. The natural frequencies will shift down due to the 
surface energy. For the first frequency, the shifting magnitude comes up to 33.5% of the original value at 
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Figure 1: The changes of first three natural frequencies with the radius of the sihca gel nanoparticle 



tq = 2.5nm, and 12.8% at rp — 5nni. Even when tq = 50nm, this ratio can also reach 1.2%. Therefore, 
the surface energy effect is not negligible in the vibration of the silica gel nanoparticle. 
On the other hand, the influence of surface energy on the higher order frequency is far smaller than that 
on the first frequency. At ro = 2.5nm, the calculation shows that the shift of third frequency caused by 
the surface energy is less than 2% of this frequency. Therefore, for the more higher frequency, the surface 
energy effect can be ignored. 




Figure 2: The changes of the first natural frequencies with the radius of the Pd nanoparticle 

Taking E = 123.6GPa, ly = 0.39, p = 12023kg/m3 and 7 = 6.0N/m [30], we investigate the radial 
vibration of the Pd particle with the surface energy effect. The results are shown in Figure 2, which 
depicts the first natural frequency of the Pd particle changing with radius. From it, it can be seen that 
the the influence of surface energy on the first natural frequency is extremely slight. The calculation shows 
that the frequency shifting caused by the surface energy is about 3% of this frequency when ro = 0.5nm. 
Therefore, the surface energy effect can be almost ignored in the vibration of the Pd nanoparticle. 



The results above indicate that, in the vibration, the influence of surface energy on the silica gel nanopar- 
ticle differs from that on the Pd nanoparticle. From the view of theory, this difference is determined by 
relative ratio the elastic modulus to the surface energy factor. Palladium is a "hard" material. Its surface 
energy factor 7 is far smaller than the elastic modulus in magnitude. This makes that the influence of 
7 on Eq. (31) is extremely slight so that it can be ignored. However, silica gel is a soft matter. At the 
nano- level, the ratio 7/Aro > 1. Therefore, the surface energy effect must be considered in the vibration 
of the silica gel nanoparticle. 

5 Summary 

A vector potential depending only on the field variables and coordinates is defined. By introducing 
the divergence of this vector potential into the Lagrangian, a Lagrangian framework is developed to 
incorporate surface energy into the theory of elasticity. Besides the Euler-Lagrange equation and natural 
boundary condition, a new boundary constitutive equation is derived from the variation of surface area 
and volume of continuum. The boundary constitutive equation characterizes the relation, on the boundary 
surface of continuum, between the Lagrangian, the vector potential and the mean curvature of surface. 
Explicit expression of the vector potential with respect to the field variable and the mean curvature is 
also determined on the boundary surface. 

Based on this framework, an elastic model with relevance to surface energy is established. In this model, 
a generalized Young-Laplace's formula is derived consistently. Compared with the existing results, the 
generalized Young-Laplace's formula depends not only on surface energy and surface curvature but also 
on the normal displacement of surface. Making use of this model, we investigate the surface energy ef- 
fect in the radial vibration of spherical nanoparticle. By calculation, some conclusions are given as follows: 

(1) At the nano-level, the influence of surface energy on the natural frequencies of a solid particle depends 
on the relative difference between elastic modulus and the surface energy factor of material in magnitude. 
The surface energy effect can be ignored in the vibration of hard nanoparticles, but for soft nanoparticles, 
it is necessary to take this effect into account. 

(2) For a soft nanoparticle, its natural frequencies will shift down due to the surface energy effect. The 
smaller the size is, the stronger the frequencies shifting. 

(3) In the vibration of soft nanoparticle, the influence of surface energy on first natural frequency is far 
greater than that on the other frequencies. For the higher order frequencies beyond the third frequency, 
the surface energy effect can be ignored. 
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